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ABSTRACT 

The compact binary system in OJ287 is modelled to contain a spinning primary black hole 
with an accretion disk and a non-spinning secondary black hole. Using Post Newtonian (PN) 
accurate equations that include 2.5PN accurate non-spinning contributions, the leading order 
general relativistic and classical spin-orbit terms, the orbit of the binary black hole in OJ287 
is calculated and as expected it depends on the spin of the primary black hole. Using the 
orbital solution, the specific times when the orbit of the secondary crosses the accretion disk 
of the primary are evaluated such that the record of observed outbursts from 1913 up to 2007 
is reproduced. The timings of the outbursts are quite sensitive to the spin value. In order to 
reproduce all the known outbursts, including a newly discovered one in 1957, the Kerr parameter 
of the primary has to be 0.28 ± 0.08. The quadrupole-moment contributions to the equations of 
motion allow us to constrain the 'no-hair' parameter to be 1.0 ± 0.3 where 0.3 is the one sigma 
error. This supports the 'black hole no-hair theorem' within the achievable precision. 

It should be possible to test the present estimate in 2015 when the next outburst is due. The 
timing of the 2015 outburst is a strong function of the spin: if the spin is 0.36 of the maximal 
value allowed in general relativity, the outburst begins in early November 2015, while the same 
event starts in the end of January 2016 if the spin is 0.2. 

Subject headings: gravitation — relativity — quasars: general — quasars: individual (OJ287) — black 
hole physics — BL Lacertae objects: individual (OJ287) 

1. Introduction elled as a binary black hole consisting of a non- 

spinning secondary black hole orbiting a more 
The BL Lacertae object OJ287 is known to have massive non-spinning primary black hole in an 

a quasiperiodic pattern of outbursts at 12 year eccent ric orbit having a perio dicity of about 12 



intervals ('SiUan paa et al. | |1988|; | Valtonen et al. | vears ()Lehto fc Valtonen Ill996l) . Further, the dou 



K)08b,). Until now, the system has been mod- pg^k structure in the light curve of OJ287 
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with the two peaks separated by 1 - 2 years, 
is interpreted as the double impact of the sec- 
ondary b lack hole on the accretio n disk of the 
primary ( Lehto &: Valtonen 19961 ). The model 



has been successful in predicting future outbursts: 
the predictions for the beginning of 1994, 1995 
and 2005 outbursts were correct within one to 
two weeks. The prediction for the 1994 outburst 
was 19 94.65 based on the average outburst in- 
terval (jSillanpaa et al. I ll988l) . The outburst had 
already begun when OJ287 was first observed dur- 
ing the fall 1994 observing season at 1994.68; 
its ar rival on time confirme d the basic period 



1996at I Valtonen fc Lehto 



icity (jSillanpaa et al 
19971 ). The probable time of the beginning of 



this outburst was at 1994.62 based on correlation 
comparison with light curves of other outbursts. 
Subsequently, an astrophysical model was con- 
structed with impacts on a rigid accretion disk 
and delays of outbursts relative to the disk cross- 
ing; the model predic ted the start of the nex t 
outburst at 1995.87 (|Lehto fc Valtonen I Il996f) . 
The observati ons gave t he best timing at 19 95.845 
(jSillanpaa et a l. 1996bl IValtonen fc Lehto] [T997). 
For the 2005 outburst, it was necessary to con- 
sider the bendi ng of the accretion disk cau sed by 
the secondary ( Sundehus et al. 19961 1997 ): when 
combined wi th the earlier model for th e radiation 
burst delay (jLehto fc Valtonen I ll996h the begin- 
ning of the 2005 outburst was expected at 2005.74. 
The best determination from observat ions gives 
2005.76 ijValtonen et al.1l2006al [2008al ). It con- 
firmed the need for relativistic precession since 
without the precession the outburst would have 
been a year later. Finally, the gravitational radi- 
ation energy loss was included i n the prediction 
for th e next outburst at 2007.70 (jValtonen 112007 . 
20081 ). Recently, the correctness of the prediction 
within an accuracy of one day was demonstrated 
( Valtonen et al. 2008bl ). The model that does not 
incorporate the gravitational radiation reaction 
effect is clearly not tenable: it predicted the out- 
burst three weeks later. The probability of the 
four major outbursts happening at the predicted 
times by chance is negligible, and the predictive 
capability of the model requires general relativis- 
tic description for the orbit of the binary black 
hole. 

We do not include the dynamical friction from 
hitting the accretion disk since the directly af- 



fected disk mass is only a few hundred solar mass, 
too small in comparison with the secondary black 
hole mass of over lO^M© to give an observable 
signal. 

In the present study, we allow the primary black 
hole to spin in the model and probe its dominant 
consequences. The additional parameters required 
to describe the dynamics of such a spinning binary 
black hole model for OJ287 can be constrained by 
recently measured historical outbursts as well as 
the more recently observed outbursts. Needless 
to say, all of these outbursts were not employed 
while constructing th e non-spinning binary black 
hole model for OJ287 ("Va ltonen eral]l2008bh . Re- 
call that the model without spin effects requires six 
well timed outbursts to provide a unique orbit and 
the usually employed outbursts are those which 
occurred in the years 1913, 1947, 1973, 1983, 1984 
and 2005. The resulting model also explains 1995 
and 2007 outbursts as mentioned above and there 
is no conflict with the observations at any other 
expected outburst times which have not been ob- 
served with adequate precision. Recently, a ma- 
jor outburst that occurre d in 1957 was identi- 
fied in the historical recor d ( Valtonen et al ] |2006bt 
Ram padara th et al.ll20071 ). It turns out that this 



identification allows one to constrain the spinning 
binary black hole model for OJ287. The next 
expected outburst in 2015 is also spin sensitive 
and may provide an improved measurement of the 
spin. 

In the literature we have found several mea- 
surements of black hole spins. Genzel et al. 
(2003) estimated the spin of the Milky Way cen- 
tral black hole as 0.52 zfc 0.2 times the max imum 
possible. Brenneman and Reynolds I ( 2006t ) esti- 
mated a spin parameter of 0.989^Q gQ2 fo^' the 
black hole in the active galaxy MCG-6-30-15 using 
relativistically-broadened spectral features from 
the inner accretion disk. There are other es- 



timates for stellar mass bla c k holes (|Liu et al. 
20081: iMcClintock et all l2006l: IShafee et al.l 120061: 



Gou et al.l 120091 ) with resulting estimates varying 



between 0.65 and 1. 

We begin by reporting the observations from 
the historical records and then explain how we de- 
scribe the PN-accurate binary orbit. Finally, we 
report the results from the timing model with dif- 
ferent spin values and present our conclusions. 
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2. Observations 

The Harvard plate stacks contains about 500,000 
glass photographic plates and as such it is the 
largest historical plate collection in the world, rep- 
resenti ng about 15% of all archival photographic 
plates (|Hudecl l2006'). The plates in the Harvard 
stacks were exposed between years 1885 and 1993 
(with a gap in 1953 - 1968; so called Menzel gap) 
and therefore they provide a possibility to gather 
optical data as old as 120 years. Plates taken 
with the same telescope belong to one plate se- 
ries. There are several such plate series at Harvard 
plate stacks. With the use of several such series, 
two of the authors (R.H. & M.B.) have managed 
to obtain more than 500 newly measured histor- 
ical data points for the blazar OJ287 and about 
3,000 data points for other 8 blazars during their 
stay at Harvard plate stacks in January-February, 
2007. In this paper we present the results of the 
measurements on 9 plates of the AC series, mea- 
surements covering the outburst of OJ287 in the 
year 1913. After conversion of our estimates to the 
standard B band, we obtain 9 new data points to 
the historical light curve of the 1913 outburst. The 
estimates of the brightness of the new data points 
were done with the naked eye equipped with a 
magnifying glass and with the help of nearby con- 
stant comparison stars of known magnitudes. We 
call this method the modified Argelander method. 
In Table [T] we give the comparison stars used for 
the estimates of the data points of the 1913 out- 
burst. Ou r comparison stars are stars 1, 2 and 7 
of Craine (jCraine Ill9"77h (hereafter C77). 



Table 1: Comparison stars used for the estimates 
of the brightness of OJ 287 around the year 1913. 
The B and V magnitudes adopted in this paper 
are the magnitudes given in C77. 
ID RA DEC B 



1 08 

2 08 

3 08 



54:29.5 
55:15.6 
55:19.0 



20:10:56 
20:09:44 
20:00:54 



12.19 
13.49 
14.13 



V 
11.52 
12.80 
13.39 



The Sonneberg Observatory has the second 
largest archive with over 300,000 photographic 
plates covering eight d ecade s beginning in 1924 



(iKroll fc Neugebauer I Il993t iKroll fc Berthold 



20051 ). The plates were taken by the Sonneberg 
Sky and Field Patrols. The Field Patrol moni- 



tored 80 fields along or near the northern Milky 
Way using astrographs and a Schmidt camera at 
limiting magnitudes up to 18. The Sky Patrol 
records the entire northern sky in two colours 
(red and blue sensitive plates) with 14 short- 
focus cameras at limi ti ng ni a gnitude of 14 - 15 



(iKroll fc Neugebauer I 1 19931 : iKroll fc Berthold 



20051 ). Scanning of the plates began in 1990. 



Thanks to Sonneberg Observatory, we obtained 
digitised photographic plates from the plate 
archives for the period 1957 - 1994. The plates 
were both blue and red sensitive, but here we de- 
scribe only results from 223 blue sensitive plates, 
as there exist more reliable comparison star se- 
quences in blue. The plates were scanned usin g 
the HP Scanjet 7400c (iKroll fc Berthold 1 12005^ . 
T he general data re duction method is described 
bv llnnis et al.l ( 20041 ) and the Windows based pro- 
gram Maxim DL was used to obtain magnitudes. 
The magnitudes were tran sformed to standard 



B sy stem as described by (jKroll fc Neugebauer 



19931 ). Thirteen nonvariable stars within 5 degrees 



of OJ287 were measured for comparison purposes. 
They were used to obtain for each plate the trans- 
formation equation from instrumental magnitude 
to B magnitude. Using the comparison stars, the 
rms error was found to vary from 0.12 to 0.20 
magnitudes from plate to plate. Some of this 
is suspected to arise from intrinsic variability of 
three of the comparison stars. Altogether we esti- 
mate the rms error of OJ287 measurements as 0.15 
magnitudes in B. The B magnitude s were trans- 
form ed to V by using B - V = 0.45 (jTakalo et al 



19941 ). Further, the V-magnitudes were converted 
to the energy flux by F = 13.18 x io(-0-4(v-i3.6)) 



where F is in mJy ([Bessell I ll979h . 

In the interval 1957 - 1978 we found 5 out- 
burst events exceeding 29 mJy at 1957, 1959 
1971, 1972 and 1973. All except the 1957 out- 
burst a re previously kno\vn and described else- 
where ( Lehto fc Valtonein Il996l : Valtonen et al 



2006br ). We have altogether 16 measurements for 



this season, out of which 3 are new (for older mea- 
surements, see Valtonen et al. 2006b) 

In Table [2] we present the newly measured data 
of the light curve of OJ287. The estimates of error 
are based on the quality of the plate. The data are 
also plotted in Figs [T] and [2] as weekly averages. 
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Table 2: New measurements (* = bad 



yr 
1912.94 
1912.94 
1912.96 
1913.00 
1913.04 
1913.08 
1913.09 
1913.20 
1913.9 
1957.149 
1957.172 
1957.180 



B 
14.00 
13.60 
14.00 
13.65 
13.20 
13.15 
13.45 
13.15 
14.15 
13.31 
13.22 
13.18 



V 
13.55 
13.15 
13.55 
13.20 
12.75 
12.70 
13.00 
12.70 
13.70 
12.86 
12.88 
12.73 



error est. 

0.3 

0.5 

0.3 

0.5 

0.3 

0.3 

0.3 
0.4 * 
0.4 * 
0.15 
0.15 
0.15 



plate) 
mjy 
13.80 
19.95 
13.80 
19.05 
28.84 
30.19 
22.90 
30.19 
12.02 
26.06 
25.58 
29.37 




Fig. 1. — The observation of the brightness of 
OJ287 in the 1912/3 season: weekly averages. The 
starting time of the outburst is estimated to be 
1912.98 ±0.02 



3. PN-accurate orbital description present 
in the model 

We invoke the 2.5PN-accurate orbital dynam- 
ics that includes the leading order general rela- 
tivistic, classical spin-orbit and radiation reaction 
effects for describi ng the tempor al evolution of a 
binary black hole (jKidder In the present 

model, where only the primary black hole is spin- 
ning, the 2.5PN-accurate inspiral dynamics of the 
binary black hole is fully provided by the follow- 
ing differential equations that describe the relative 
acceleration of the binary and the precessional mo- 
tion for the spin of the primary black hole. The 
2.5PN-accurate equations of motion can be writ- 




1956.95 1957 1957.05 1957.1 1957.15 1957.2 1957.25 1957.3 
Year 

Fig. 2. — The observation of the brightness of 
OJ287inthe 1956/7 season: weekly averages. The 
starting time of the outburst is estimated to be 
1957.08 ±0.03 



ten schematically as 



X 



fx 



= Xo+XipN +XSO+XQ 
+ X2PN + X2.5PN , 



(1) 



where x = Xi — X2 stands for the center-of-mass 
relative separation vector between the black holes 
with masses mi and m2 and Xq represents the 
Newtonian acceleration given by Xq = — -^^^a;; 
m = uii + 7712 and r = \x\. The PN contributions 
occurring at the conservative IPN, 2PN and the 
reactive 2.5PN orders, denoted by asipjv, X2pn 
and X2.5PN respectively, are non-spin by nature. 
The explicit expressions for these contributions, 
suitable for describing the binary black hole dy- 
namics, were derived for the first time in the har- 
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monic gauge (jDamour 1119821 ). These are given by 



reads 



XiPN 



X2PN 



Gm 



-2(2 + r/) 



Gm 



(2) 



+ (1 + 3t])v^ - ^T]P] - 2(2 - ?;)ft)| , 



Gm 



n 



(12 + 297?)(G2»)^ 

+??(3 - 4r/)u4 + ^r;(l - 3r])r^ 
r]{3 - Ari)v^P - ^t]{13 - 47y)(^) 

Ar])v'^ - (4 + 41?7 
-3ry(3 + 2r])P 



(3) 



-(2 
[^(15 4 



X2.5PN 



_ _8_ Gjmfi) 



15 



[9„2 _,_ 17^] ^„ (4) 



[3i 



iGml 



where the vectors n and v are defined to be 
n = x/r and v = dx/dt, respectively, while 
r = dr/dt = n ■ v, v = \v\ and the symmet- 
ric mass ratio rj — mi m2/m^. The leading order 
spin-o rbit contributions to x, app earing at 1.5PN 
order ([Barker fc O'Connell iflQTSh . reads 



Ji. _ Gm /Gm\ ( 1 

xso - — [—) y- 



4 



(5) 



12 [si ■ (n X v)] 



[9 + 3 Vl-4?7) 



(n X si 



7 + VI -477 



{v X Si 



where the Kerr parameter x ^^'^ the unit vector 
Si define the spin of the primary black hole by 
the relation Si = Gm\x^i/c and x is allowed 
to take values between and 1 in general rela- 
tivity. Further, the above expression for xso hn- 
plies that the covariant spin supplementary condi- 
tion is employed to define the center-of-mass world 
line of the spinning co mpact object i n the underly- 



ing PN computation ( Kidder 1995f) . Finally, the 
quadrupole-monopole interaction term xp^ enter- 



ing at the 2PN order (jBarker k O'Connell IITgTSh . 



Xc 



2 3G^ mini 



5(n • Si)2 - 1 



n 



-2{n ■ Si)si 



(6) 



where the parameter q, whose value is 1 in general 
relativity, is in troduced to test the black hole 'no 
hair' theorem j|Thornelll980t IWex and Kopeikin 
Il999t |Wil]ll2Q08|). The precessional motion for the 



primary black hole spin is dominated by the lead- 
ing order general relativistic spin-orbit coupling 
and the relevant equation reads 



dsi 

f Gmri\ /7 + ^Jl — Arj 



\2c^ r2 



i + yr^T^ 



(n X d) . (7) 



It should be noted that the precessional equation 
for the unit spin vector Si enters the binary dy- 
namics at IPN order, while the spin contribution 
enters x at the 1.5PN order. 

The main consequence of including the leading 
order spin-orbit interactions to the dynamics of 
a binary black hole is that it forces both the bi- 
nary orbit and the orbital plane to precess. More- 
over, the orbital angular momentum vector, char- 
acterising the orbital plane, precesses around the 
spin of the primary in such a way that the an- 
gle between the orbital plane and the spin vector 
Si remains almost constant (roughly within ±0.°5 
in our model). This is illustrated in the Figure 
[31 The spin-vector itself precesses drawing a cone 
with an opening angle of about 12 degrees (see 
Figure HI). 

The precessional period for both the orbital 
plane and the spin of the binary, provided by | f2 1 , 
is about 2400 years, while employing the ahead 



inferr ed orbital parameters for OJ287 (Valtonen 



ay 
in I 



20071 ) . The relevant information required to evalu- 
ate |r2| are the primary mass mi = 1.8 x Mq, 
the secondary mass m2 = 1.2 x 10^ Mq, the or- 
bital period of 9 yr (without the redshift factor), 
and the eccentricity at apocenter ~ 0.66. Fur- 
ther, the orbital inclination relative to the plane 
of symmetry of the spinning black hole ( and pre- 
sumably relative to the accretion disk of the pri- 
mary as well) is taken to be 90 degrees. The man- 
ner in which the orbital plane precesses around Si 
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20000 




-20000 I ' ' ' ' ' ' ' 1 

-20000 -15000 -10000 -5000 5000 10000 15000 20000 

X 

Fig. 3. — Precession of the orbital plane due to the 
post-Newtonian spin-terms in a model in which 
the spin vector is kept constant (for clarity of the 
figure). The system is seen from the direction of 
the spin axis which in our model is in the orbital 
plane and this situation persists with the precision 
of approximately half a degree when the plane it- 
self precesses. The time interval illustrated is 750 
years during which the plane precesses about 180°. 
At the same time the periastron precesses much 
more rapidly and thus we see the orbit alternately 
from the direction of aphelion, perihelion or side- 
ways, which causes the apparent different lengths 
of the small BH excursion from the big BH (which 
is located at the origin of the figure). 

is displayed in Figure [31 In the next section, we 
probe its consequences to the disk crossing tim- 
ings. Contrary to the illustration in Figure [3l the 
spin precession is taken into account in the timing 
experiments. 

4. Timing experiments 

In previous work (e.g. Valtonen 2007) the tim- 
ing experiments were performed using the 1913, 
1947, 1973, 1983, 1984 and 2005 outbursts as fixed 
points. They determine 5 time intervals whereby 
it is possible to determine uniquely 5 parame- 
ters. They do not include the secondary mass, 
and neither can we determine the Kerr parame- 
ter X and the no-hair parameter q without con- 
sidering more fixed points. Here we also make 
use of the 1995 and 2007 outbursts which have 
been observed with high time resolution (see Fig- 
ure [5]) and they are therefore suitable for further 
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Fig. 4. — The circulation of the ascending node 
( fl) of the orbit (straight line) in the coordinate 
system in which the initial spin vector si points to 
the direction of the z-axis. The quantities plotted 
are fi/lOO in degrees and the symbols Sx and Sy 
correspond to the x and y components of the spin- 
direction vector Si in degrees. The plot for thus 
illustrates (partly) the same behaviour as FiglH] 

refinement of the model. In addition, the newly 
discovered 1957 outburst is taken as another fixed 
point. The observed outburst times are listed in 
Tabic [3] Note that we have used the beginning of 
the outburst as defined by the light curve rather 
than a linear extrapolation from the steepest ris- 
ing slope. This change makes the outburst times 
a little earlier than was previously assumed; e.g. 
the 1983 outburst timing is now 1982.964 rather 
than 1983.00 (see Figure [5]). With this full set of 
outbursts no solutions were found unless the pa- 
rameter X is in the range 0.2 - 0.36. The majority 
of solutions cluster around x — 0.29, with a one 
standard deviation of about 0.04. 

Figure [6] shows the distribution of points, each 
representing a solution, in four panels. Each panel 
displays the outburst time as a function of x- Pan- 
els (a), (c) and (d) are essentially scatter dia- 
grams, demonstrating that the solutions cover the 
allowed range rather well. On the contrary, panel 
(b) shows a strong correlation between the time 
of the 2015 outburst and x- The range of possi- 
ble outburst times extends from early November 
if X = 0.36 to late January 2016 if x = 0.2. 

The timing experiments give a unique solution 
for the system parameters: precession in the or- 
bital plane per period A0, masses mi and m2 of 
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Beginning of 1995 outburst 
at 1995.842 



2.5 3 3.5 4 4.5 5 5.5 
Day in November 1995 




10 10.5 11 11.5 12 12.5 13 
Day in September 2007 



16.5 

16 - 
15.5 

15 - 
14.5 

14 



Beginning of 1973 outburst 
at 1972.94 



8 10 12 14 16 
Day in December 1972 



Beginning of 1983 outburst 
at 1982.9645 , 



Day in December 1982 



0.3 - I 

0.25 -• 

0.2 -* 

0.15 -' 

0.1 - 
0.05 



1.06 
1.04 
1.02 
1 

0.98 
0.96 
0.94 
0.92 



0.2 0.24 0.28 0.32 0.36 



0.2 0.24 0.28 0.32 0.36 



0.46 1— 

0.44 - 

0.42 - 

0.4 -• 

0.38 - 

0.36 - 
0.34 



0.2 0.24 0.28 0.32 0.36 



0.36 
0.34 
0.32 

0.3 
0.28 
0.26 
0.24 
0.22 

0.2 
0.18 
0.16 



0.2 0.24 0.28 0.32 0.36 



Fig. 5. — The beginning of the best timed out- 
bursts at 1973, 1983, 1995 and 2007. For 1973 
and 1983 the points are daily averages, for 1995 
averages over 8 hours, and single observations for 
2007. The curve is a model fit for a spherical ho- 
mogeneous bubble which becomes instantaneously 
transparent at the initial moment, given inside the 
figure. 

the two black holes, spin parameter initial phase 
(j)o, initial apocenter eccentricity eo, 'no- hair' pa- 
rameter q and time-delay parameter td- The re- 
sults are shown in Table 21 

Here the precession rate is defined as the aver- 
age change of the apocenter phase angle over the 
last 150 yrs. The time delay parameter depends on 
the structure of the accretion disk, and obtains the 



value of unity in the model of iLehto &: Valtonen 
(|l996l ). Its value is related to the thickness of 



the acc retion disk, as explained in I Valtonen et al 
(|2006bl ). 



Because we are using 8 time intervals to fit 8 pa- 
rameters, the fit has no degrees of freedom. The 
tolerance limits of Table 4 are a consequence of 
having a tolerance in the fixed points of the out- 
burst times (Table 3). We have carried out fur- 
ther experiments trying to evaluate how well de- 
termined these tolerance limits are. In case of q 
values outside the range of Table 4 were found in 
a limited part of the initial value space. Thus the 
error limits for this quantity are likely to be one 
sigma error limits. For other variables we were 
not able to get values outside the range except 
by changing the astrophysical model, as reported 



Fig. 6. — Variation of the expected outburst tim- 
ing as a function of black hole spin for four differ- 
ent outbursts. Results from the various acceptable 
solutions are represented by points. 

Table 3: Outburst times with estimated uncertain- 
ties. These are starting times of the outbursts. 

1912.980 ± 0.020 

1947.283 ± 0.002 

1957.080 ± 0.030 

1972.945 ± 0.012 

1982.964 ± 0.0006 

1984.130 ± 0.005 

1995.842 ± 0.0015 

2005.745 ± 0.012 

2007.692 ± 0.0015 



below. 

As part of the solution we obtain the list of 
all past and future outburst times with their 
uncertainties. For example, the well recog- 
nized outbursts in 1959, 1971 and 1994 are 
timed at 1959.213 ± 0.002, 1971.1265 ± 0.002 and 
1994.6085 ± 0.005, respectively. In all these cases 
data are missing at the crucial time of expected 
rapid flux rise, and thus these predictions cannot 
be verified at present. 

The value for the primary mass is within errors 



the sa me as was obtained by iLehto fc Valtonen 
(Il996l) . mi = 1.71 ± 0.15 • IO^OM^. Also the sec- 



TOi = i.(iznu.itj - lu iviQ. 

ondary mass is surprising ly close to the value of 



Lehto fc Valtonen I (199 



TO2 



1.36 • 10** M, 



0: 



when adjusted for the 'current' Hubble constant 
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Table 4: Solution parameters. 





39.1 ±0.1 




(1.83 ±0.01) • lOi°M0 




(1.4 ±0.1) • 1O*^M0 


X 


0.28 ±0.08 


00 


56P.5 ±1.°2 


eo 


0.6584 ±0.001 


q 


1.0 ±0.3 


td 


0.75 ±0.04 



of 72 km/s/Mpc and for the 5.6 mJy outburst 
streng th, as observed in 2007 (jValtonen et al. 
[2008b). This mass value is based on the as- 
trophysics of the impact and the strength of 
the maximum signal. The previous values of 



Lehto fc Valtonen j (|l996l ) were 80 km/s/Mpc and 



5.1 mJy which correspond to the secondary black 
hole mass 1.0 • 10^ A/0. The updated astrophys- 
ically determined secondary mass is within the 
error limits of the dynamically determined mass. 
Considering that the two ways of measuring the 
mass are completely independent, the agreement 
is a good indicator of the correctness of the basic 
model. 

We may ask what is the likelihood of finding 
a quasar like OJ287 in an all sky search for peri- 
odic quasars. One may approach this issue from 
a theoretical point of view in cosmological evolu- 
tion models, or as an observational problem. In 
the former case it has been calculated that many 
binary black holes are expected in quasars, e.g. 
Volonteri et al. (2009) estimate that there should 
be about 10^ sub-parsec binaries up to redshift 
z = 0.7 among quasars with absolute z-magnitude 
less than -22, and mass ratio not exceeding 100 
over the whole sky. For a system like OJ287 there 
is something like 25% chance to be found at a red- 
shift below 1 in cosmological simulations (private 
communication by an anonymous referee). 

When considering the observational situation, 
we note that at its faintest OJ287 has been ms = 
18; since the jet of OJ287 is strongly beamed to- 
ward us, we use this base magnitude in compar- 
ing with other quasars. There are about 2-10'* 
quasars brighter than this limit in the sky (Arp 
1981). Many of these quasars may host binary 
black holes (Comerford et al. 2009). Thus poten- 
tially there are about lO'* bright (at OJ287 base 
level) binary quasars in the sky. The rate of bina- 



ries in quasars is somewhat uncertain at present. 

Let us again consider the special case of OJ287. 
Vestergaard & Kelly (private communication) esti- 
mate that there is about a unit probability of find- 
ing a broad-line quasar brighter than i-magnitude 
18 where the black hole mass exceeds lO^'^Af0 in 
a search covering the whole sky. The probability 
of finding one hundred such quasars is about 25%. 
This is an extrapolation from the study of SDSS 
(DR3) quasars (Vestergaard et al. 2008), com- 
plemented by the BQS sample. The effect of the 
bright limit {i — 15) in the SDSS catalog which 
reduces the number of observed nearby massive 
quasar black holes, has been corrected for. Alto- 
gether, the estimate based on observations is not 
inconsistent with the estimate from the cosmolog- 
ical evolution models, considering the large uncer- 
tainties in both. 

According to our model, OJ287 is in the stage 
of final inspiral before the merger of the two black 
holes. Its remaining lifetime of about 10^ yr is 
only a fraction 10^'^ of the expected time of black 
hole binary evolution of about 10^ yr (Volonteri et 
al. 2009). On the other hand, OJ287 is strongly 
beamed toward us. Even though the likelihood 
that we happen to be looking straight at the beam 
of a quasar is only 7~^/4 where 7 is the Lorenz 
factor of the beam, the intensity of the source is 
magnified by 7'^, and the probability of detection 
in a fiux limited sample goes as the 1.5-power of 
the intensity. The joint probability that OJ287 en- 
ters the samples of both optically bright and radio 
loud quasars (the selection criteria for the Tuorla 
monitoring program) is raised by several orders 
of magnitude above the corresponding probabil- 
ity for unbeamed quasars. At its brightest, OJ287 
is about 5 magnitudes brighter than at its mini- 
mum light. This increase in detection probability 
counters the decrease in detection probability due 
to the limited lifetime of its currect stage of evo- 
lution. Based on these considerations, it is not 
surprising that a short period binary quasar with 
the primary heavier than has been dis- 

covered. 

There is nothing special about other parameters 
of OJ287 either. The large mass ratio is necessary 
for the stability of the accretion disk. The moder- 
ately high eccentricity is what is expected at the 
OJ287 stage of inspiral of binary black holes of 
large mass ratio (Baumgardt et al. 2006, Matsub- 
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ayashi et al. 2007, Sesana et al. 2008, Iwasawa 
et al. 2009). The value for td implies a '--^ 0.1 in 
the standard disk model, an acceptable value. In 
the next section we will discuss the main targets 
of this paper, the values of x and q. 

We have performed numerical experiments to 
study the stability of the above listed solution. 
In the first instance, we added the conservative, 
non-spinning 3PN corrections to x describing the 
relativistic binary black hole orbital dynamics 
([Mora fc Willll2004l ). The timing experiments em- 
ploying the improved binary black hole dynamics 
lead to orbital solutions that are essentially close 
to what is listed in Table 3. The only notica- 
ble difference is the roughly 1.5% increase in the 
estimated mass of the primary black hole. This 
result suggests that the employed PN-dynamics, 
relevant for our binary black hole model, is in 
the convergent regime. The increase in the esti- 
mated mass of the primary can be qualitatively 
explained by evaluating the period of periastron 
precession Tproc at 2PN and 3PN orders with th e 
help of Eqs. (20) in 'Mem mesheimer et all ( 2004 ). 
Using a non-spinning binary black hole model for 
OJ287, we obtain Tpicc ~ 87 years at the 2PN 
level and an increase in the primary mass by 1.5% 
reduces it to ~ 86 years. The similar estimate in- 
volving the non-spinning binary black hole model 
at the 3PN order gives Tprec ^ 83 years indicat- 
ing that our timing experiments are qualitatively 
consistent with our rough theoretical estimates 
for Tproc- The quantitative differences may be 
attributed to the fact that in the timing exper- 
iments, we employed spinning binary black hole 
model that decays under gravitational radiation 
reaction. We would like to point out that it is the 
data set spanning roughly 100 years that makes 
Tpicc a good diagnostic tool in analyzing our tim- 
ing experiments. Let us also note that at present 
the 3.5PN reactive contributions to x do not in- 
fluence our estimates in any siginificant way. This 
is because its contribution to the rate of change of 
the binary black hole orbital period is ~ 0.00014 
and this is within the error limits of our analysis. 

The second test involves varying the astrophys- 
ical parameters of the model. The free parame- 
ter td already tests the different acceretion disk 
structures. Any errors in the model for the evolu- 
tion of the radiating bubble affect primarily this 
parameter, not others. There are two more astro- 



physical parameters which have fixed values in our 
model; here we will briefly discuss their influence 
on the solution. The disk level above or below 
the mean level is affected by the tides raised by 
the secondary. The effect of tides on the impact 
times was calculated by Valtonen (2007), and it 
is described by an analytical function of the im- 
pact distance. We have tested possible errors in 
this formula by varying the disk tidal level by a 
constant factor. It turns out that the level can 
vary anywhere between 40% and 200% of the stan- 
dard value without affecting our result. No so- 
lutions were found if the disk level correction is 
outside this range. Another fixed parameter in 
the model is related to the asymmetry of the im- 
pact with respect to the central plane of the disk. 
First the secondary impacts the upper surface of 
the disk, initiates the evolution of the radiating 
bubble, and then one disk crossing time later it 
impacts the lower surface of the disk and initi- 
ates the evolution of the second radiating bubble 
(Ivanov et al. 1998). Since we observe the disk 
from one side only, we see one or the other of the 
two bubbles. Thus our model must include the 
delay by one crossing time when we observe the 
bubble from the 'lower' level. We derive the value 
1.05 ± 0.15 crossing times for the delay from the 
illustrations of Ivanov et al. (1998). Since in the 
double impacts the radial speed changes its sign 
relative to us between the approaching part of the 
orbit and the receding part of the orbit, every sec- 
ond outburst must include the one crossing time 
delay. What happens if the delay is not exactly 
one crossing time? We have varied the delay by 
a constant factor, and found that if the delay is 
within the range of 0.85 — 1.20 crossing times, it 
does not affect our result. It is unlikely that the 
delay would in fact fall outside this range. How- 
ever, if the delay is for some unknown reason much 
shorter than one crossing time, solutions are found 
down to the level of 0.55 of the crossing time, but 
then the spin value is greater than in our standard 
model. At the extreme, the spin can be as high as 
X = 0.43. As the correlation of Figure 6b shows, 
this possibility can be eliminated (or verified) in 
2015. The same correlation still holds even if the 
delay is much less than unity. 
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5. Conclusions 

The binary black hole in OJ287 is modelled to 
contain a spinning primary black hole with an ac- 
cretion disk and a non-spinning secondary black 
hole. Using PN-accurate dynamics, relevant for 
such a system, we infer that the primary black hole 
should spin approximately at one quarter of the 
maximum spin rate allowed in general relativity. 
In addition, the 'no-hair theorem' of black holes 



( Misner at al.l [l973,) is supported by our model, 
although the testing is possible only to limited pre- 
cision. The solutions concentrate around q = l.O 
with one standard deviation of 0.3 units. These 
results are achieved with the help of new data on 
historical outbursts as well as using the most re- 
cent outburst light curves to gether with the ti ming 

IqqtI) 



model for OJ287 outbursts ( Valtonen 



To 

obtain the above estimate for the Kerr parameter 
X and the no-hair parameter q, we have assumed 
that the spin of the primary black hole is initially 
(150 yrs ago) aligned with respect to the spin of 
accretion disk. Further, a polar orbit is assumed 
as in previous work. The main justification for 
the polar orbit is the sudden fade of OJ287 over 
a wide range of wavelengths in 1989, presumably 
because of an eclip se of the jet by the secondary 
( Takalo et al.lll990( ) . The polar orbit was also used 
to justify the rather modest rise of flux in 1995, 
possibly because the secondary moved parallel to 
the disk axis after the disk impa ct and blocked 



off som e of the optical emission (jValtonen et al 
2006bh . 



What happens if these conditions are not sat- 
isfied? Sundehus et al. (1997) carried out sim- 
ulations with a number of different inclinations 
between the disk and the orbit, and found that 
the inclination makes no difference. This is not 
surprising since the impacts occur along the line 
of nodes, and the timing of the impacts is not a 
function of the impact angle. In this work we have 
tried varying the spin angle relative to the disk, 
and found that within reasonable limits (we con- 
sidered inclinations up to 10°) the results do not 
change. Again this is what is expected since the 
line of nodes circulates slowly in comparison with 
the orbital angular motion. 

The reason why we would expect at least ap- 
proximate alignment between the black hole spin 
and the disk spin is the Bardeen - Peterson effect 



which tends to align the two spins. However, the 
alignment time scale is much slower (about 10^ yr, 
Lodato & Pringle 2006) than the spin precession 
time (Fig. 4), and thus the disk does not follow 
the black hole spin exactly but only on average. 

In this scenario, we have a unique solution and 
also a unique prediction for the next OJ287 out- 
burst, expected in 2015. We should then be able 
to judge the correctness of the present solution. 
Note that an outburst is not expected in 2015 in 
any simple extrapolation from past observations, 
as it is well before the the average 12 yr cycle is 
due, and thus it is a sensitive test of the general 
model as well as a test for the spin of the primary 
black hole. 

In the meantime, there is an additional observa- 
tion which supports the spin value that has been 
derived here. OJ287 h as a basic 46±3 day peri- 
odicity ( Wu et al. 2006[ ). which may be related to 
the innermost stable orbit in the accretion disk. 
However, since we presumably observe the accre- 
tion disk almost face on, and there is an m = 2 
mode wave disturbance in the disk, this is likely 
to refer to one half of the period. Considering also 
the redshift of the system, and the primary mass 
value given in Table 4 this corresponds to th e spin 
of X = 0.35 ± 0.06 (|McClintock et al.ll2006^ . The 
uncertainty of 0.06 units is related to the width 
of the trough in the structure function of the flux 
variations. 

It will be interesting to explore the consequence 
of allowing the secondary black hole to spin. This 
will require including the leading order spin-spin 
interactions, appeari ng at the 2P N order, to the 



inspiral dynamics ( Kidder I I1995I) . We are also 



planning to probe the effect the higher order con- 
tributions to the spin-orbit coupling, entering the 
binary black hole bina ry dynamics at the 2.5PN 
order (jPave et al.ll2006f ). 

There are at least 5 additional outbursts in the 
historical record which have not yet been been ac- 
curately observed. If new data of these outbursts 
are found it will open up the possibility of im- 
proving the model and studying higher order ef- 
fects. Searches of plate archives are planned for 
this purpose. 
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